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O ' Abstract 



o ' The elastic hadronic amplitude is calculated using the nonperturbative light-cone 

p. 



Mh. dipole representation for gluon bremsstrahlung. The data for large mass diffraction 



demand a two-scale structure of light hadrons: the gluon clouds of the valence quarks 
I^ ' with the size of ~ 0.3 fm and the hadronic size ~ 1 fm. The presence of the two 

j^ ■ scales unavoidably leads to a specific form for the total hadronic cross section which 

consists of a steeply rising oc s (A = 0.17 it 0.01) term related to gluon radiation, 
and a large constant term originating from soft interactions which does not induce any 
gluon emission. Our calculations reproduce well the total cross sections and elastic 
slopes Q. To further test the model, we analyze the elastic pp and pp differential cross 
sections and extract the partial amplitudes in the impact parameter representation. 
The Pomeron trajectory as a function of the impact parameter is only slightly above 
one for central collisions, but steeply grows towards the periphery. The model predicts 
correctly the shape and energy dependence of the partial amplitude at all impact 
parameters. 
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1 Introduction 

1.1 Rising total cross sections 

The increase of the total hadronic cross section at high energies is well known since the 
ISR experiments in the early 70's. This discovery came just in the time when the Regge 
theory had its conjecture. The simple idea to shift the intercept of the Pomeron pole above 
one, ap(0) = 1 + A, leads to contradiction with the unitarity restrictions, in particular, 
the Froissart bound is violated. It took special efforts to formulate a self-consistent Regge 
scheme |0] in which unitarity is restored via Regge cuts and without violation of energy 
conservation (warnings of which had been given in [Q] ) . 

The assumption that the Pomeron which governs the hadronic elastic amplitude at high 
energies is a Regge pole has no theoretical justification beyond simplicity. It faces problems 
interpreting data from HERA which demonstrate that A substantially increases with Q^ 
(see, however, [Q, P). 

The ensemble of data on hadronic elastic scattering at small t, i.e. total and differential 
cross sections, slopes and ratios of real to imaginary part of the forward elastic amplitudes 
can successfully be fitted by many phenomenological models based on varieties of quite 
different assumptions aimed to fit the data (see, for example |^). Even the simple parame- 
terization s^ with exponential i-dependence for the elastic amplitude describes well the data 
at small t [0, §]. However, the imaginary part of the forward elastic amplitude is connected 
by the unitarity relation to the total inelastic cross section, in other words, the Pomeron is 
the shadow of inelastic processes. Unavoidably, one should study the dynamics of inelastic 
collisions to understand the forward elastic scattering 0, rather than guessing its analytic 
form which is only mildly restricted by general principles. 

The total cross section for a highly virtual photon interacting with a proton measured 
in deep-inelastic lepton scattering (DIS) can be estimated using perturbative QCD if the 
photon virtuality Q^ ^ Qq {Qq ~ IGeV) and the energy s ^ Sq (sq ~ IGeV"^), but 



X = Q^/s <^ 1. Depending on the approximations used, two models for the hard Pomeron 
are known: the BFKL |^ and the double- leading- log DGLAP (see in |jlO|, |lll). The rising 



energy dependence is interpreted perturbatively as caused by gluon bremsstrahlung with 
growing phase space for radiated gluons. The total cross section is predicted to rise steeply 
with energy as is confirmed by data from HERA. The energy dependence, parameterized as 
g^eftiQ ) reveals the exponent to increase with Q^ up to Ag// ~ 0.5. 

In terms of the QCD light-cone dipole approach one can treat DIS at small x as an 
interaction of a tiny size, ~ 1/Q, quark-antiquark fluctuations surrounded by a gluon cloud 
which is much larger (logarithmically) than the qq pair. 

1.2 Soft interaction limit: the two scales for light hadrons 

A new scheme for performing explicit calculations for the interaction of light hadrons has 
been suggested in M. It exploits the smallness of the gluon correlation radius which has 



been estimated in many approaches. In particular, the model developed in ref . |T^ , extends 
the perturbative methods of the light-cone QCD to the nonperturbative region introducing 
the light-cone potential into the Schrodinger equation for the Green function that describes 
the propagation of a quark-gluon interacting pair. The interaction potential fixed by the 
data for large mass soft diffraction turns out to be rather strong leading also to a short 
separation tq = 0.3 fm between gluons and the source (a quark or a gluon). This result is 
confirmed by the recent analysis [J13l of HERA data for diffraction which leads to an even 
smaller estimate tq ~ 0.2 fm (but with large uncertainties). Thus, a proton looks in the 
infinite momentum frame like three valence quarks surrounded by small gluon clouds as 
illustrated in Fig. [^. 



Such a two-scale structure of light hadrons appears not only in the model [0 treating the 
vacuum fluctuations as Weizsaker- Williams gluons. The smallness of the gluon clouds of the 
valence quarks is confirmed by the study of the gluon formfactor of the proton employing 



QCD sum rules |T^. The Q^ dependence of the formfactor turns out to be rather weak 



Figure 1: A skeleton of three valence quarks in the proton, surrounded 
by gluon clouds of much smaller size than the mean quark separation. 

corresponding to a small radius of the gluon distribution which was estimated at the same 
value tq ^ 0.3 /m. The small gluon correlation radius ~ 0.3 /m appears also from lattice 
calculations |T^. It is also predicted by the liquid instanton model |[T^, |T^] and is related 



to the instanton size po ~ 0-3 frn. The experimental observation of a small cross section 



for large mass soft diffractive dissociation has led to a small value of tq |T2|, |1^] and this, 
quite in general, can be taken as the confirmation of the small size cloud of any kind of 
gluonic vacuum fluctuations surrounding the valence quarks. These are usually referred to 
as constituent quarks although nothing specific (about mass, additivity, etc.) is assumed 
beyond the simple statement that the clouds of vacuum fluctuations dressing the valence 
quarks are much smaller than the mean hadronic radius. 

Of course the transverse size of the gluonic spots increases with energy since the weights 
of higher Fock components grow as powers of ln{s) [|l], [T^. Such a behavior is specific 
of gluonic fluctuations. Nevertheless, the mean size of the fluctuation clouds is still small 
compared to the radii of light hadrons in the energy range of modern accelerators (see [Q 
and below). The ratio of the constituent quark radius tq to the mean interquark separation 
Rh squared serves as a small number. Correspondingly, one should single out two different 
contributions to the total inelastic cross sections. 

The first one is due to the soft interaction which is unable to resolve the structure of 
the constituent quarks and excite them. This contribution can be treated as the cross 



section of interacting hadrons made of structureless valence quarks which can be viewed 
as the skeleton of the hadrons. One can try to evaluate it using either the naive two- 



gluon exchange approximation |T8|, IT^, BO, ETl O, |22[ or more sophisticated nonperturbative 



approaches like string crossing and rearrangement ||2^, or the interaction of overlapping 
Wilson loops in the stochastic vacuum model [^, |25[. This part of the cross section o"o(i?/i) 



is controlled by the mean interquark separation Rh and is independent of energy since the 
size of the quark skeleton of the hadron is constant. 

The second contribution to the total cross section comes from the semi-hard interaction 
able to resolve the small size of the constituent quarks and to excite them giving origin 
to gluon radiation. This cross section is proportional to the size of the constituent quark, 
oc Tq. The radiation of each new gluon leads as usual to an extra power of ln{s) which 
exponentiates to a cri(ro)s'^ dependence. However, the energy independent term cri(ro) 
needed for the exponentiation of these logs is rather small (oc Tq) and can not match the 
large term ao{Rh)- 

Thus, we arrive at the following general structure of the total cross section which corre- 
sponds to a two-scales scheme for the hadronic structure, 

CTtoi = c^o + 0-1 (ro) ( — ) , (1) 

\Sq/ 

where do = ao{Rh) — o"i(ro). Parametrically, ao{Rh) ^ cri(^o)- 

The double scale structure of light hadrons, (i.e. small constituent quarks versus large 
interquark separation), leads to the structure of the total cross section Eq. (|l|), rather more 
complex than the usually assumed overall behavior oc s^. Of course, in the spirit of the 
leading-log (s) approximation, one can neglect the constant term ctq as s — *> oo, but then A 
should not be compared with the experimental data available in the energy range where ctq 
gives an important contribution. The effective slope of the energy dependence may, in fact, 
be substantially smaller, 

Ae// = (l - ;p^) A . (2) 

In fact, in [0 it is found that Ae// ~ A/2. 



1.3 Outline of the paper 

This paper is organized as follows. In section we calculate the cross section of gluon 
radiation in high energy hadronic interactions. It is shown that valence quarks contribute 
additively to the radiation cross section due to the short-range correlation of the radiated 
gluons. Moreover, the sum of the multigluon radiation cross sections which depend on energy 
as powers of Ins exponentiates in the leading-log approximation to the energy dependent 
power s^. The exponent proportional to the running QCD coupling a^ turns out to be 
rather large, A = 0.17 ± 0.01 compared to what is believed to be demanded by the present 
data for total cross section. However, gluon radiation contributes with a rather small factor 
proportional to Tq ?^ 1 mb. This energy dependent fraction of the total cross section is fully 
predicted. The large energy independent part of the cross section is due to the interaction 
of the valence quarks with no gluon radiation. This cross section is related to the large 
hadronic size, rather than to tq, and cannot be evaluated perturbatively. Although it can 
be estimated in models, e.g. as it is done in the stochastic vacuum model [^, ^, the 



uncertainty of such calculations is too large, and we prefer to treat do as a free parameter, 
which turns out to be the only unknown of the model. 

The rising total cross section eventually violates the Froissart-Marten bound at very 
high energies, but the partial elastic amplitude at small impact parameters is already very 
close to the limit imposed by unitarity and may easily break it down. The procedure of 
unitarization of the elastic partial amplitudes is described in section ^ We use the standard 
quasi-eikonal model, but we compare also with a different QCD motivated approach. 

The model is analyzed with respect to the total cross section data in section ^. The only 
parameter of the model, (Jq, can be fixed by comparison with the data at any chosen energy. 
Then, the energy dependence is predicted in a good agreement with the data. The slope of 
forward elastic scattering needs no new parameters and is also well predicted. 

In the standard Regge phenomenology the energy dependence of the total cross sections 
and of the elastic slopes are controlled by the intercept and slope of the Pomeron trajectory. 



respectively, which are independent parameters. However, one may expect them to be 
correlated since in QCD the cross sections depend on the hadronic sizes due to color screening 
|I8|, pJ|, pn, p6[. An attempt to incorporate this property was made recently in p^. Here we 



develop this approach treating more consistently the phase space for the radiated gluons. 

The comparison with the data turns out to be most effective in the impact parameter 
representation. First of all, the radius of interaction exposes explicitly in this case. Secondly, 
unitarity imposes severe restrictions on the elastic partial amplitude for central hadronic 



collisions |28|, which slows down the energy dependence of the partial amplitude [29, 2S]. 
Thirdly, the color dipole representation in QCD introduced in [^ became a popular tool 
to study high energy QCD dynamics in DIS, Drell-Yan reaction etc., since color dipoles are 
the eigenstates of the interaction at high energies. In this respect, the impact parameter 
representation is suitable for a direct comparison of the data with a dynamical model (see for 



instance [0). And, last but not least, the shape of the amplitude in the impact parameter 
space is related to the shape of the amplitude as function of momentum transfer in a wide 
range of t, rather than only in the forward direction. 

In section ^ we analyze the available high energy data from ISR and SppS for pp and pp 
elastic scattering to extract the partial elastic amplitude in the impact parameter represen- 
tation. We follow the procedure suggested by Amaldi and Schubert p9[ who performed a 



similar analysis of ISR data and concluded that the total cross section rises due to peripheral 
interactions while the partial amplitude for central collisions is energy independent. This is 
usually treated as a manifestation of unitarity saturation. However, the parameterization 
of the amplitude used in |2^ was based on the geometrical scaling model which assumes 
that the ratio of the total cross section to the elastic slope is independent of energy. This 
assumption unavoidably leads to a constant partial amplitude at zero impact parameter (see 
Fig. i. 

The geometrical scaling is known to be broken beyond the ISR energy range, therefore 
we perform the analysis differently, in a less model dependent way. We fit the t-dependence 
of the cross section independently at each energy assuming no correlation between different 



energies, except the normalization which is adjusted to fit the energy dependent total cross 
section and ratio of real to imaginary parts of the forward elastic amplitude. The t-dependent 
imaginary part of the amplitude arising from the fit is then Fourier transformed to the impact 
parameter representation at each energy. The 6-dependence of the partial amplitudes found 
this way is very close to what our model predicts. Not only the shape of the partial amplitude 
is well reproduced, but also its development as function of energy. 

In section ^ we compare the data for the partial amplitudes at different energies and 
conclude that they hardly vary at 6 ^ 0, but rise steeply with energy at large 6 > 1 fm. 
The effective Pomeron trajectory is a steeply rising function of the impact parameter. Our 
model correctly predicts this dependence. 

The results of the paper are summarized in section |^. Further evidences of the large 
value of A suggested by data on diffraction in DIS and particle production at mid rapidities 
in soft hadronic collisions are reviewed. 

2 Excitation of valence quarks: nonperturbative gluon 
radiation 

To calculate the energy dependent total cross section one should sum up the various con- 
tributions of different Fock components of the incoming hadron. To avoid double-counting, 
we sum the cross sections o"„ of physical process of radiation of n gluons. 



O". 



hN _ sr^ _hN 
tot 

n 



T.^'n''- (3) 



The lowest Fock component (n = 0) of a hadron contains only valence quarks. For the 
sake of simplicity we assume the beam hadron to be a meson, the generalization to a nucleon 
is simple and is done below. 

The contribution to the total cross section corresponding to the interaction without any 



gluon radiation has the form, 



cr, 



hN 



doq I d^R 



^iia^R) ag{R) 



(4) 



Here the valence quark wave function of the hadron, \l/|g(ag, R), depends on the transverse 
q — q separation R (see Fig. |^) and the fraction a^ of the hght-cone momentum of the pair 




R-r 



Figure 2: A cartoon for the Fock state \qqG) in the impact parameter 
plane. 



carried by the quark. The energy independent Born cross section of interaction of a large 
qq dipole with a nucleon cr^^^R) cannot be calculated perturbatively since the separation R 
is large. It can not be adjusted directly to the experimental data since the data include the 

3n i 
we treat a^^^ as a free parameter. 



contribution from gluon bremsstrahlung leading to the energy dependence of a^ . Instead 



The next contribution to the a^^ comes from radiating a single gluon. The radiation 
is possible only due to the difference between the amplitudes for the qq and qqG Fock 
components, otherwise the interaction does not alter the combination of Fock states and 
they remain coherent, i.e. nothing new is produced. Another way to explain this is to say 
that the interaction with the target can free the gluon fluctuation only if it resolves it, i.e. 
discriminates between the interaction amplitudes for two Fock components \qq) and \qqG). 



The contribution to the total cross section corresponding to radiating a single gluon 



reads pl §, 



a 



hN 



fdagfd^R\^'^giR,ag)f J -^ J d^r 



aG<l 



X 



^g-G(i? + r, ac) agiR + f ) + %G(r, ac 



'<(-) 



- Re v^^clr, ac) ^.ciR + r, ac) <(i? + r) + <(r) - <(i?) 



,Af/ 



.A^/ 



(5) 



Here ac is the fraction of the hadron momentum carried by the gluon which is assumed 
to be small; the notations for the radii are obvious from Fig. |^. The first and the second 
terms in the curly brackets correspond to the emission of the gluon from the quark and the 
antiquark respectively, and the third term to the interference between them. 

The nonperturbative wave function for a quark-gluon Fock component was derived in 



| 12| . Neglecting the quark mass the wave function reads, 

2i [du e* ■ r 



'^qGir^ac] 



exp 



2rl 



(6) 



a<l TT V 3 r" 

where e is the polarization vector of the massless gluon. The mean separation tq = 0.3 fm is 
related to the nonperturbative light-cone potential describing the quark gluon interaction. 
It is fixed by the data on large mass diffractive dissociation corresponding to the triple- 
Pomeron limit. The mean quark - gluon separation r^ is much smaller than the distance R 
between the quarks. Therefore, one of the qG wave functions in (|^) can be neglected leaving 
a factor 2 (a factor 3 in the case of A^A^ scattering) since both q and q can radiate the gluon. 

At small separations, r ~ tq, the dipole cross section cr^g{r) can be evaluated pertur- 
batively and the approximation cr^(r) = C r"^ can be used. The two-gluon approximation 
gives for the factor C ~ 2.3 using an effective gluon mass rriG = 0.15 GeV to incorporate 
confinement, and a^ = 0.4 (see below). 

Thus, the contribution of the \qqG) Fock component to the total cross section summed 
over polarization of the radiated gluon takes the form. 



CTi 



hN 



4a f s\ 

In — 

37r \soJ 
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Cr' 



(7) 



Here ln[s/so] = ln[(Q;G)max/(«G)mm] originates from the integration over ac in @, where 
{c(G)min = 2/srQ ^ {lGeV'^)/s, but (Q;G)maa: IS iH defined. It should be small enough, say 
~ 0.1 to make sure that the quark-gluon wave functions in (^ is independent of ac- Then, 
assuming that the quark carries a fraction one third of the proton momentum, we estimate 
the value Sq ~ 30 GeV, which we use in what follows. 

The radiation of each new, n-th, gluon can be treated as the radiation by an effective 
quark, which is the valence quark surrounded by n — 1 gluons. It should be resolved by 
the soft interaction with the target as being different from the radiation oi n — 1 gluons. 
Therefore, it provides the same mean cross section 9Cro/4 and the factor Aag/Sn as in 
(0). This is illustrated in Fig. ^ in the 1/Nc approximation, i.e. replacing each gluon by 
a qq pair. According to the general prescription \^2\ , |12| of the light-cone approach, the 




q^^^r 8 

Figure 3: Radiation of the second gluon in the leading-log (s) approxi- 
mation as seen in 1/Nc approximation when each gluon is replaced by a 
qq pair. Solid quark lines correspond to the final state, dashed - lines 
correspond to the initial state configurations. 



radiation cross section in the impact parameter representation is proportional to the total 
cross section of a colorless system made of all the final state partons (solid lines in Fig. |^) 
plus the initial state partons replaced by antipartons (dashed lines). Since the radiation 
of gluons with ac ^ 1 does not affect the impact parameter of the radiating quark, all 
the solid and the corresponding dashed lines have the same impact parameters. Therefore, 
those quark-antiquark (solid-dashed) pairs which are color neutral do not contribute to the 
cross section. The only quark contributing is the one that radiates the last gluon and that 
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changes color in a color-octet state with the antiquark as shown in Fig. ^. Thus, the total 
cross section of the multiquark configuration is reduced to one for the octet-octet dipole 
with mean separation tq, i.e. dCr^/A. 

In conclusion, the n-th term in (Q) reads (for a single valence quark). 






4ao , / s 



<^s 



In 



-I n 



Icrl. (8) 



-^'o' 


+ 3 


4 ° 


(-) 




A = 


4 a, 






Stt 




~< 


9 

-^0 4^ 


^2 



Summing up the powers of logarithms in (^) we arrive at the following expression for the 
total cross section, 

cyZ = ~<' + ^-.Crl (-) , (9) 

where 

Am 

(10) 

(11) 

Since each of three valence quarks can radiate (see (^) the second term in (^ acquires 
factor 3. 

The structure of Eq. (|^) refiects the physical input illustrated in Fig. |l|, as discussed in 
the Introduction (section 1.2). The energy dependence of the cross section is related to the 
excitation of the small spots (constituent quarks) inside the hadron, while a large energy 
independent contribution corresponds to the soft interaction of the valence quarks skeleton 
of the hadrons leading to no gluon radiation. These two parts of the cross section cannot 
match to provide a power dependence, s^, as a common factor, as it is usually assumed to 
the the case in the so-called soft Pomeron approach ||^, |21], |2^ . This is the origin of the more 



complicated form we find i.e. of the structure of eq. (^) or, which is the same, of eq. (|1]). 
The energy dependent part of the cross section, i.e. the second term in (^, is suppressed 
by the smallness of Tq and is expected to be relatively small at medium-high energies, but 
grows more steeply with energy than the overall (7^4. Note that such a structure of the total 



cross section was also suggested in [33|, however with a different physical motivation. It was 



found to fit well the data on total cross section (see below). 
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The power A in (j^) is related to the Pomeron intercept, ap(0) = 1 + A, and can be 
predicted using ( [Tol) provided the QCD couphng a^ at virtuahty ~ l/rg is known. In 
Gribov's theory of confinement ||3^, ^ the radius of a constituent quark is at the borderhne 
between the perturbative and the nonperturbative QCD regimes. At larger distances, the 
chiral symmetry breaks down and pseudo-Golstone pions emerge^. At smaller distances. 



perturbative QCD is at work. The corresponding critical value of as is |34, 35 




a, = — \l-J-\^0A3. (12) 



Another way to evaluate a.^ is to average the running QCD coupling weighted with the 
transverse momentum distribution of radiated gluons, 

, ^ Jo aK^as\i^T) dilnaG)dk-:^ . . 

JO "'^T d{lnaG)dk'^ 

where the transverse momentum distribution of gluons radiated in quark-nucleon interaction 



is given by Eq. (130) of p 



The standard phenomenological way to extend ^^(/cr) down to small values fcy — > is 
to make a shift in the argument, k^ ^ k^ + k^. The value k^ ^ 0.25 GeV^ was estimated 
3^ using dispersion techniques |^ of higher twist effects in hard reactions. 



m 



We evaluated ([13|) using the dipole cross section <Jqq{p) oc 1 — exp{— p^ / p^) which is 

proportional to p^ at small p, but levels off at large separations. The nonperturbative 

quark-gluon interaction taken into account in [|T^] is very important, since it squeezes the 

quark-gluon fluctuations down to a mean size tt ~ tq substantially increasing the mean 

transverse momenta of radiated gluons. Correspondingly, the mean value {as) turns out to 

be rather small. For the parameter po varying within a reasonable interval 0.3 < po < 1 fm 

we found the mean coupling varying between (a^) = 0.38 — 0.43 which agrees well with the 

critical value Eq. (p!2D. Substituting the central value of (a^) into ([T0|) and using the interval 
^Thc Pomeron properties at large distances related to pion loops have been employed in the model 
suggested in M 
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as the uncertainty for a^ we get, 

A = 0.17 ±0.01. (14) 

This value is about twice as large as the value 0.08 usually believed to be required by 
data. This value, however, can not be compared directly with the soft Pomeron intercept 
where it is assumed that the whole cross section is proportional to s^. The value of ( |T^ ) 
is just the second term of (|^) and the overall energy dependence is much less steep as a 
consequence of the large value of the constant term ctq^. We will show later (see Section ^) 
that the predicted energy dependence of the cross section (^ is in a good accord with the 
data and corresponds to an effective value A^ff ~ 0.1. First of all, however, we should take 
care of unitarity since the cross section (H) violates the Froissart bound at large s. 

3 Impact parameter representation, unitarization 

Although the total cross section stays well below the Froissart-Martin bound up to the 
present highest energies, the partial amplitude at small impact parameters demonstrates 
a precocious onset of the unitarity restrictions which are already important in the energy 
range of existing accelerators. 

The imaginary part of the partial amplitude corresponding to the total cross section (|^) 
can be decomposed into the terms related by unitarity to the radiation of different number 
n of gluons, 

Im7p(s, 6) = ^Im7„(s, 6) , (15) 

n 

where 7„(s,6) is the partial elastic amplitude which depends on the energy and on the 
impact parameter b. Upon performing the Fourier transform of the t-dependent elastic 
amplitude, the integral over b gives the corresponding term a^^ in (H). We assume that the 
t-dependence of the lowest Fock component which is related to the spatial distribution of 
valence quarks is given by the product of the electromagnetic formfactors of the colliding 
hadrons (we confine our considerations to pp and pp collisions) F^{t). For simplicity we use 
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the standard dipole form for the proton fornifactor Fp{t) = (1 — t (r^^)/12)~^, where (r^^) 
is the mean charge radius squared related to the slope of elastic scattering of the valence 
quark skeleton hj Bq = 2 (r^^)/3. 

We keep the same t-dependence for higher Fock components in (|^) corresponding to 
gluon radiation by a projectile valence quark interacting with the target proton; the slope 
of these components, however, should increase linearly with the number of radiated gluons 
due to their random walk in the impact parameter plane with a step ~ Tq for the radiation 
of every new gluon, 

Bn = I (rl,) + ^ . (16) 

The Fourier transform of the square of the dipole formfactor leads to the following shape 
for the partial amplitudes ||26|| , 



(T'i^(s) 3 



lmYnnb,s) = ^^y'K,iy), (17) 

where y'^ = {4b'^/Bn)^, K^{y) is the third order modified Bessel function and a^^ are given 
by (^) and (§). The normalization of the partial amplitude is fixed by the relation, atot = 
2 Jd%lm-f{b,s). 

The partial amplitude ([15|) rises with energy and eventually would lead to a violation 
of the unitarity bound Im7(s, 6) < 1, unless unitarity corrections are introduced. Unfor- 
tunately, this is not a well defined procedure since different recipes can be found in the 
literature. 

The simplest known way to restore unitarity is to eikonalize the partial amplitude (p!5|). 



Imrp(6, s) = 1 — exp — Im7p(6, s 



'U 



At very high s this amplitude approaches the black disk limit 0, Imrp(s, b) -^ 0[i?^(s)— 6^], 
with radius, R{s) = ro Aln(s/so). Correspondingly, at high energies 

Aln(f)»^ (19) 
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all hadronic cross sections reach the maximal universal energy growth allowed by Froissart- 
Martin's bound, 

at^{s)^2nrlAHn'[£l . (20) 

The eikonalization procedure ([l8| ) would be suitable if the the incoming hadrons were 
eigenstates of the interaction [Q. Hadrons, however, are subject to diffractive off-diagonal 
excitation, and the eikonal form of unitarization should be corrected in a way similar to Gri- 



bov's inelastic corrections [^ for hadron-nucleus cross sections. The lowest order unitarity 
correction in (|I^) comes from the quadratic term in the exponent expansion of r{b,s). It 
has to be modified using the AGK cutting rules to include single diffraction, 



1 2 

ImTp = Im7p - - (lm7p) [l + D{s)] + 0{-fl) 



(21) 



where D{s) = o"sd(s)/(Je;(s) is approximately 0.25 in the ISR energy range and decreases 
slightly with energy oc 5"°°^ [g|, |0]. Indeed, aeii-s) = al^{s)/{l6TTBei) oc s°•^ but the 
energy dependence of the diffractive cross section is rather flat (due to stronger unitarity 
corrections), cTsdis) oc s°'°^. Asymptotically, as s ^ oo, D{s) vanishes since aei{s) oc In^s 
and Csdis) oc Ins. 

The inelastic corrections to higher order terms in the expansion (|T8|) are poorly known. 
A simple way to keep (^) and to include diffraction into the higher terms is to modify (|T^) 
as, 



Imrp(6, s) 



1 



1 + D(s] 



1 - cxp -(l + D{s)^lm'jp{b,s) I 



(22) 



which is known as quasi-eikonal model [^ 



A more consistent way of unitarization suggested in [Q prescribes to use the eikonal 
expression (jT^) in terms of the color dipole cross sections and then average it over the 
transverse separations of all partons. Unfortunately, this procedure is simple only if the 
dipole cross sections depend quadratically on the separation parameter, which is definitely 
incorrect for cr^(-R) in (^). If it were true one would have -D(s) = 1 which would exceed four 
times the experimental value. We have tried this unitarization prescription as well and found 
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that the results still agree with data pretty well. Nevertheless, we use for further applications 



Eq. ( P2|) since it explicitly exploits experimental information, correctly reproduces the lowest 



unitarity correction (pT]), and is rather accurate within the energy range of interest. 



4 Comparison with forward scattering data 



All the parameters in Eqs. (|^, (^) and ( P^ are known, except the Born cross section 
o"o. Although it may be estimated in various models, none of these is sufficiently reliable 
and we choose to determine it so as to best reproduce the data. As soon as the absolute 
normalization of the total cross section is fitted at some energy and the parameters a^ is 
fixed, the energy dependence can be predicted. For this comparison we selected the data 
| 42| for crfo( at ^/s = 546 GeV as being the most precise. In addition, this energy value is 



high enough that we neglect the Reggeon contribution. We calculate a^ot from 

atot = 2 Jd%lmT{b,s) , (23) 



using (p2D and fix our only unknown parameter at o"o = 39.7 mb. 

Now we are in the position to predict the energy dependence for the total cross sections 
and compare it with pp and pp data. The result is depicted by the dashed curve in Fig. |^ is 
in a good agreement at high energies, but somewhat off the data at medium high energies. 
This is not surprising since the (secondary) Reggeon contribution is still missing and this is 
well known to be important at medium high energies. 

In order to improve the description of the data one should, therefore add the contribution 
of leading Reggeons with intercept aji{0) ~ 1/2. This should be done directly in the partial 
elastic amplitude, 

Im T{s, b) = Im Tp{s, b) + Tr{s, b) [l - Im Tp{s, b)] . (24) 

The Reggeon term is suppressed by the absorptive corrections which have the same origin 
as those which slow down the energy dependence of the diffraction cross section mentioned 
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Figure 4: Data for total pp (closed circles) and pp) (open circles) cross 
sections P^/ at y/s > lOGeV. The dashed curve shows the predicted 
energy dependence for the net Pomeron contribution whose normaliza- 
tion is fixed by the y/s = 546 GeV data Ji^/ . The solid curves, (bottom 
for pp and upper for pp), represent the results corrected for Reggeon 
contribution which is fitted to the data. 



above. The Reggeon term is parameterized as, 



lmTR{s,b) 



(^R 



^tiBr{s 



Sq 



ajj.(0)-l 



exp 



2Br{s 



(25) 



and Br = R\ + 2a^lns. We fixed the standard values the parameters ai?(0) = 0.5 and 
a'^ = 0.9 GeV^^, but fitted Rj^ = 3GeV^^. We also fitted the normalization factors which 
are very different for pp and pp because of the (approximate) exchange degeneracy. We 



found a^ = 17. 8 mb and a^ = 32.8 mb. The result of the fit, shown in Fig. ^ by solid 
curves demonstrates a good agreement with data. As anticipated, the Reggeon corrections 
are important only in the ISR energy range and below. 

Since o"o and all the Reggeon parameters are now fixed, we can predict the slopes of elastic 
pp and pp scattering (both the absolute values and energy dependence). We calculate the 
slope using the relation, 

Beiis) = Ub^) = — fd%bHmT{b,s) . (26) 

Once again, the results shown in Fig. ^ in comparison with data for pp and pp scattering 
demonstrate a good agreement. Although we had some freedom in the choice of the proton 
formfactor and of the proton charge radius, this affects only the absolute value of the slope. 
The energy dependence is fully predicted. Since it describes the data well, we correctly 
predict the effective Pomeron slope a'^^^ ~ 0.25 GeV^^. 

Note that, often, phenomenological fits treat total cross sections and slopes as controlled 
by different parameters. In these cases, one cannot predict the energy dependent slope even 
if the total cross section is known. 

The radiation of every new gluon leads to an expansion of the gluon cloud by a "step" 
(5(r^) K, (0.3 fmf'. Eventually, the initial approximation of a small gluon cloud inside a 
large hadron will break down. This, however, will happen only at very high energies. The 
mean number of gluons in a quark (n) = Aln(s/so) is quite small, (n) = 0.5 — 0.8 at ISR, 
about {n) ^ 1.5 at SppS and reaches (n) ?^ 2 at the Tevatron. Therefore the mean radius 
of a constituent quark is still rather small and our approximation remains quite valid and 
we should expect it to break down only at very high energies which are well beyond the 
range of present accelerators. Of course the radius of an constituent quark, i.e. the radius 
of the gluon cloud depends also on a reference frame, e.g. in the cm. it is twice as small as 
in the rest frame of the target. 
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Figure 5: Data for the elastic slope ki4l and our predictions. The upper 
and bottom curves and, correspondingly, the open and full circles belong 
to pp and pp, respectively. 



5 Elastic scattering data analyzed in the impact pa- 
rameter 



The partial amplitude ( [2^ ) has nontrivial s- and 6-dependences. It is nearly energy inde- 
pendent for central collisions, but steeply grows with energy on the periphery as it was first 



found in the analysis of the data by Amaldi and Schubert p9| . These properties are averaged 
out and hidden in the total or differential elastic cross sections. To extract information about 
the shape of the partial elastic amplitude in the impact parameter representation from the 
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data on elastic pp and pp scattering, we follow the procedure used in |29|. However, to make 
the analysis less model dependent we fit differential elastic cross section data independently 
at each energy, thus, no model for energy dependence is involved. The geometrical scaling 



model used in |^ assumes that the total cross section is proportional to the slope of the elas- 
tic differential cross section, 0-^4(5) oc Bpp{s). This relation is, however, a result due to the 
accidental closeness of the Regge model parameters characterizing the energy dependence 
of the cross section, atoti^) oc s'^ ~ 1 + A Ins and of the slope, B^^^s) = Bq^ + 2aplns, 
where a'p ~ 0.25 GeV~^ and Bq^ ^ 7.5 GeV""^. Indeed, the effective Pomeron intercept 
A = ap(0) — 1 ^ 0.08 is close to the ratio 2a'p/ Bq^ ^ 0.067. Obviously, geometrical scaling 
may occur only in a restricted energy range (namely, in the ISR energy range used in ||29[| ) 
and it had been predicted |@] to break down at higher energies. This was confirmed later 
by the SppS and Tevatron data. 

Since we are interested in the Pomeron part of the elastic amplitude, the data for differ- 
ential cross section of elastic pp and pp scattering selected for the analysis cover the wide 



range of high energies including /S'-RQ, and SppS ||5]. We do not include the data from 



the Tevatron since they are available only in too narrow a range of t which is not sufficient 
for Fourier transformation. We parameterize the imaginary and real parts of the elastic 
scattering amplitude in a model independent way as, 

Im/(t) = ^a,e^^*; (27) 

Re/(t) = ^Qe'^»*, (28) 

where Oj, 6,, q, di are parameters to be fitted. The amplitudes are related to the cross 
sections as, 

[Re/(t)]' + [Im/(t)]^ (29) 



dt 

atoi = 4v^Im/(0) . (30) 



^We rely upon the relative normalizations of differential cross sections measured in different f-intervals 
at the same energy given in p9(| 
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To make the normalization of the differential cross section data more reliable, first we 
perform a common fit of the pp and pp total cross sections with the same Pomeron part 
as function of energy. Then we adjust the normalizations of the differential elastic cross 
section data to the optical points, i. e. demand that 4 y/n J^o^i = o'tot at each energy. The 
data 1^ for p{s) ratio of the real to the imaginary parts of the amplitudes at t = were 



also involved in the analysis. We fit these data by a smooth energy dependence and demand 
then J2ci = P J2(^i for each energy included in the analyses of differential cross sections. 
We performed two variants of fit, 

• variant I: both Ci and C2 are used as free parameters in the fit; 

• variant II: C2 = in (p^). 

The data in the fit and the results in variant I are depicted in Fig. R. 



As soon as the parameters in (|27|) and (^) are found we can calculate the partial 



amplitude in the impact parameter representation at each energy , 

r(&) = ^/Ae^^^V(-g'), (31) 

where g is the transverse component of the momentum transfer, t ~ — g^. It is normalized 
according to (|23|). 

A few examples of our results for Imr(6) corresponding to variant I (full points) and 
variant 11 (open points) are shown in Fig. |^ with spacing 0.2 fm in the impact parameter 
and for a few energies. The errors are calculated using the error matrix resulting from the 
fit. 

One can see that at 6 = the amplitude nearly saturates the unitarity limit and hardly 
changes with energy, while at larger impact parameters the amplitude grows quite substan- 
tially. 

Our predictions including the Pomeron contribution and Reggeon part are compared 
with the data in Fig. |^. The Reggeons shown by dashed curves are calculated for pp and pp 
interactions for ISR and SppS data respectively. Their contribution is quite a small fraction 
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Figure 6: The differential cross sections of elastic pp and pp scattering 
at different energies. The first five panels show the pp data from ISR 
j^] , the last one the pp data p^/ from SppS. The curves show our fit 
Eqs. (]2^j-(|2^j in variant I. 
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of the full partial amplitude represented by the solid curves. The agreement between the 
data and our predictions is remarkably good, especially if we recall that the Pomeron part 
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Figure 7: The imaginary part of the partial amphtude Im r(6) as function 
of the impact parameter at different energies. The first three panels 
correspond to the ISR, the last one to the SppS data. The curves show 
our theoretical prediction with Eq. ( ]2^ ) using the parameters given in 
Table 1 and obtained fitting the t = data (total cross sections and 
slopes). 
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has no free parameters, except one, do, adjusted to the total cross section measured at one 
energy ^/s = 546 GeV [Q. Both the predicted shape of the partial amplitude and its energy 
development are confirmed by the data. 



6 Pomeron trajectory in the impact parameter space 

The partial elastic amplitude rises with energy faster for peripheral than for central collisions. 
The energy dependence of Imr(6, s) at different values of the impact parameter is shown in 
Fig. ^ for variant I. One can see by the eye that the upper curves corresponding to central 

imr(Vs) 
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Figure 8: Imr(6) plotted at various values ofb = 0.0, 0.2, 0.4, ... 2.8 fm 
as function of energy. The values and error bars correspond to Fig. [^. 
The lines correspond to the fit with power dependence on energy. 
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collisions are nearly horizontal, while the bottom ones representing peripheral collisions rise 
steeply with energy. The curves show the results of the fit to the data for energy dependence 
of the partial amplitude at different values of the impact parameter by the expression 

Im r(6, s) = To s'^(^) + Im Tnib, s) , (32) 



in which the Reggeon term ( P5| ) is calculated with the parameters fixed by the fit to total 
cross section data. The exponent A (6) varies with the impact parameter and is fitted to the 
data for what concerns its energy dependence in each bin of b as it is shown in Fig. ^ We 
ignored the data at ^/s = 62 GeV since they are too much off the smooth interpolation of 
the data at lower and higher energies. The results of the fit for A (6) are plotted in Fig. ^ 
as a series of black points corresponding to each of the lines in Fig. ^. Open points show 
the results of the fit corresponding to variant II. The error bars are determined by the error 
matrix of the fit. Note that these values of A (6) correspond to the Pomeron contribution 
since the Reggeon part is sorted out. 

Although A(6) is very small ~ 0.03 for central collisions, {i.e near h = 0), it increases 
dramatically (nearly by one order of magnitude) at large b. Thus, the data show that the 
energy dependence of the total cross section originates mainly from peripheral interactions. 
This confirms the observation of |2^ . 



The systematic uncertainty of our analysis is related mainly to the choice of parameter- 
ization for the elastic amplitude. The difference between variants I and II can be treated as 
a characteristic estimate of this uncertainty. There is no significant difference between the 
two solutions. 

Our predictions plotted as solid curves in Fig. |^ can also be translated into values for 
the effective exponent A (6), 



din 



rp{s,b) 



using the theoretical amplitude ( P^ with the same parameters already determined. The 
results are shown as a solid curve in Fig. 0. It agrees well with the data. 
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A(b)=d ln(r)/d ln(s) 
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Figure 9: The exponent A (6) found by the fit to each point of Fig. ^ 
with power dependence on energy at each value ofb. The black and open 
points correspond to the fits with paranieterizations I and II respectively. 
Our predictions with Eq. ( |3^ ) are shown by solid curve. The dashed 
curve demonstrates prediction of a single Regge pole model without any 
unitarity corrections. 



The dashed curve shows the prediction of the simplest parameterization ||^, H for the 
elastic amplitude with power s - and exponential t - dependences for both the Pomeron 
and Reggeon terms. Although the unitarity corrections are neglected, this parameterization 
is indeed quite successful in describing total cross sections and elastic slopes [^. One can 
see, however, that its agreement with the data in Fig. is quite poor. It overshoots the 
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data for central (lack of unitarization) and peripheral collisions and has quite a different 
6-dependent shape. Nevertheless, one should not interpret the smallness of A at 6 = as a 
manifestation of saturation of unitarity. 

7 Summary, discussion and outlook 

We present the first successful quantitative dynamical description of small angle elastic 
scattering of light hadrons. The key points of our approach are, 

• The data for diffractive gluon radiation (large mass diffraction) demand a small trans- 
verse separation tq ^ 0.3 fm between the radiated gluon and the valence quark [|T2|, ^ . 

• A new regime is found which allows explicit calculations: the gluon clouds of valence 
quarks are much smaller than the hadronic size, r^/ R\ <ti 1. This is a different limiting 
case compared to DIS where the gluon cloud is much larger than the qq separation. 

• The interference between amplitudes of gluon radiation by different quarks is sup- 
pressed as exp{—Rl/rl) leading to an additivity of the valence quarks in the part of 
the total cross section related to gluon radiation. 

• Since gluon radiation is controlled by the semi-hard scale Q ~ 2/ro it can be evaluated 
perturbatively. The radiation cross section is suppressed by the small factor r^ ^ 1 mb, 
but steeply rises with energy oc s^, where A is given by Eq. (pjlj) . 

• The part of the total cross section do related to soft collisions without excitation of 
the valence quarks is large since it is controlled by the large hadronic radius. It is 
independent of energy and may cause deviation from quark additivity. 

• The two-scale structure (ro versus Rh) of light hadrons unavoidably leads to the spe- 
cific form (|l]) of energy dependence for the total cross section. The terms with and 
without gluon radiations are governed by different scales and cannot match in order 
to exponentiate into a common factor s^ for the whole cross section. 
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• While the second, energy dependent term in Eq. (P can be evaluated perturbatively, 
none of current models can estimate reliably the first constant term ctq. We treat it 
as a free parameter and fix it by the normalization of total cross sections. Then we 
predict the energy dependence of the total cross section and the forward elastic slope 
in good accord with data. 

• As a further rigorous test of the model, we perform a model-independent analysis of 
available high-energy data for the elastic differential cross sections and extract the 
partial amplitudes in impact parameter representation. Our model describes well the 
observed b- and s— dependences of the partial amplitude. 

• We also extract the Pomeron trajectory ap{b) = 1 + A(6) in the impact parameter 
representation. A(6) is very small at 6 = as a result of unitarity saturation, but rises 
by an order of magnitude for peripheral collisions, in good accord with our predictions. 

Concluding, the strong interaction of radiated gluons is vital for present approach. It 
squeezes the gluon clouds of valence quarks and allows to apply perturbative QCD to cal- 
culation of the radiation cross section. It is worth emphasizing that it is not legitimate to 
mimic these nonperturbative effects introducing an effective gluon mass tug ~ 0.7 GeV as is 
frequently done in the literature. Indeed, only the light-cone gluons interact nonperturba- 
tively during their long lifetime. However, the t-channel Coulomb gluons cannot be treated 
on the same footing as the light-cone ones, their lifetime is always short. These gluons are 
massless and can propagate far away. To incorporate the confinement one should assign 
only a small effective mass rriG ~ ^qcd to the t-channel gluons. Making these t-channel 
gluons as heavy as those on the light-cone would suppress the factor C in (^ and the term 
o"i in ([1|) by nearly order of magnitude in contradiction with data. 

Note that our results explain the surprisingly high effective Pomeron intercept A ^ 0.2 
observed in diffractive DIS. It is known that diffraction is dominated by soft interactions 
even at high Q^ and one could expect about twice smaller value. However, large mass 
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diffraction related to diffractive gluon radiation corresponds to tlie second term in Eq. @, 
i.e. value of A given by (|14|) . 

Altliougli data for diffractive dissociation in soft liadronic collisions were used to fix the 
strength of the nonperturbative gluon interaction, one can further test the model performing 
a similar analysis of data for the differential cross section of single diffraction in the impact 
parameter representation. Such an analysis is in progress and will be published elsewhere. 

In this paper we concentrate on calculation of the elastic hadronic amplitude related 
via unitarity to inelastic processes like gluon radiation. It is natural to extend the test of 
the model comparing directly to data for multiparticle production. In particular, the AGK 
cancelation [§] of unitarity corrections leads to the inclusive cross section in central region of 
rapidities rising with energy as s^ where A has its genuine value not disturbed by unitarity 
corrections. Such an analysis of data performed in M^ has led to a surprisingly similar 



conclusions as ours. Namely, (i) the data cannot be described by the energy dependence 
oc s^ in the whole energy range, but demand an additional constant term; (ii) the fit to 
data resulted in A = 0.17, exactly what is predicted by our calculations. 



Note that the model suggested in [Q which describes multiparticle production in terms 
of energy independent string fragmentation and rising with energy mini-jet contribution also 
goes along with the basics of our model. 
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